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We review recent results on the effect of a specific type of quenched disorder on well 
known 0(m)-vector models in three dimensions: the XY model (3DXY, m — 2) and 
the Ising model (3DIS, m — 1). Evidences of changes of criticality in both systems, 
when confined in aerogel pores, are briefly referenced. The 3DXY model represents 
the universality class to which the A-transition of bulk superfluid 4 He belongs. Ex- 
periments report interesting changes of critical exponents for this transition, when 
superfluid 4 He is confined in aerogels. Numerical evidence has also been presented 
that the 3DXY model, confined in aerogel-like structures, exhibits critical exponents 
different from those of bulk, in agreement with experiments. Both results seem to con- 
tradict Harris criterion: being the specific heat exponent negative for the pure system 
(q3dxy — —0.011 < 0), changes must be explained in terms of the extended crite- 
rion due to Weinrib and Halperin, which requires disorder to be long-range correlated 
(LRC) at all scales. In numerical works, aerogels are simulated by the diffusion limited 
cluster-cluster aggregation (DLCA) algorithm, known to mimic the geometric features 
of aerogels. These objects, real or simulated, are fractal through some decades only, 
and present crossovers to homogeneous regimes at finite scales, so the violation to 
Harris criterion persists. The apparent violation has been explained in terms of hid- 
den LRC subsets within aerogels [Phys. Rev. Lett., 2003, 90, 170602]. On the other 
hand, experiments on the liquid-vapor (LV) transition of 4 He and N2 confined in aero- 
gels, also showed changes in critical-point exponents. Being the LV critical-point in 
the O(l) universality class, criticality may be affected by both, short-range correlated 
(SRC) and LRC subsets of disorder. Simulations of the 3DIS in DLCA aerogels can 
corroborate experimental results. Experiments and simulations both suggest a shift in 
critical exponents to values closer to the SRC instead of those of the LRC fixed point. 
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1. Harris criterion in brief: original and extended 

Whether the presence of disordered impurities affects the critical behavior of an ideal 
system, or not, has been the task of numerous works through years. Since Harris' seminal 
work [ m , a robust theoretical background has emerged to establish conditions for the 
relevance of disorder to phase transitions, which concern the criticality of the original 
pure system as well as the geometrical features of the disordered distribution of defects [ 
El El HI IS] On the other hand, many experimental contributions in the last two decades 
posed interesting questions about criteria of relevance, supplying results that challenge 
predictions made by previous theoretical works. Amongst all, those on the superfluid 
transition of 4 He in light aerogels are rather intriguing [ H3 |S1 H3 El . 

Changes of critical exponents for the so called A-transition, when superfluid 4 He is 
confined in aerogel pores, were reported repeatedly from late 80's through late 90's, with 
almost the same question left open: do these results violate Harris and/or other relevance 
criteria? The question arises because aerogels are homogeneous (non-correlated) beyond a 
finite scale [El El, and given that the specific heat exponent (a) is negative for this tran- 
sition, exponents should not change after Harris criterion. Many numerical and theoretical 
works have since emerged to explain these contradictory results [ rjj El El El El El • 

According to Harris criterion [^, short-range correlated (SRC) disorder is irrelevant for 
the critical behaviour of any d-dimensional pure system which undergoes a second order 
phase transition with a correlation length exponent ^ pure > 2/d. The criterion was shown 
valid if disorder presents a correlation function 5(r) and, after Josephson hyp er scaling [E]) 
otpure < 0. This original work was improved several years later by Weinrib and Halperin 
(WH)[[2], who established a more general criterion of relevance: even if a < for the pure 
system, it will change critical exponents if disorder is "correlated enough" . A proper defini- 
tion of this long-range correlated (LRC) disorder makes use of the impurity-impurity corre- 
lation function g{r) = (n(r)n(O)) and its long-range scaling exponent (g(r) ~ r~ a , r — > oo). 
Depending on how fast the tail of g(r) decays, the extended criterion reads as follows: 



2 — du pure > for a > d (1-1) 
2 — aup Ure > for a < d. (1.2) 

In both cases disorder is relevant. Case 
(jl.lj) extends the particular definition 
of short-range correlated (SRC) disorder, 
made by Harris, to disorder distributions 
with a "fast-decaying" tail (a > d). The 
second case (jl.2|) is valid even if v pure > 
2/d. Based on one— loop expansions in e = 
4 - d < 1 and 5 = 4 - a < 1, WH also ar- 
gued that, the correlation length exponent 
would be u Lnc = 2/a, independent on the 
internal dimension (m) of the order param- 
eter. Prudnikov et al. [E] performed more 
accurate estimates to deduce m-dependent 
corrections, i.e., under the same kind of 
LRC distribution of defects (same a), their 
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Figure 1. Correlation length exponent for 0(m)- 
vector models (m = 1,2,3) with LRC defects, 
depicted as functions of df = d — a/2. WH 
m-independent result plotted for comparison 
(dashed line). Data from Table V in Ref. [Ej- 
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two- loop field theoretical expansions give different critical exponents for different 0(m) 
systems. Fig. ^ resumes graphically what they obtain for the correlation length exponent 
^lro a t disorder correlation exponents 2 < a < 3 and m = 1, 2, 3, quite different from the 
m-independent plot predicted by WH 

2. Phase transitions in aerogels: experiments 

In this section, we review some important experimental results about critical systems 
confined in aerogels. In 1988, Moses Chan et al. [ 6] reported on the influence of quenched 
disorder on the superfluid (SF) transition of confined 4 He, using three types of porous glass: 
Vycor, aerogels and xerogels. They measured the temperature dependence of the relative 
superfluid density, which scales as p s /p ~ & near the critical point, being t = |T — T c \/T c 
the reduced temperature, and the exponent £ = v due to hyper universality [03^5- In the 
first case ( He— Vycor), they observed no change in the critical exponent from £ ~ 0.6705 
of bulk SF 4 He, a fact clearly explained in terms of Harris criterion, provided that the pure 
critical system exhibits ci negative specific heat exponent {oipure — 

-0.0105) [US]. Internal 

microchannel structures of Vycor, randomly oriented and randomly distributed, present 
scattering intensities that exhibit a peak near r max ~ 20nm, and fall off exponentially 
above this peak, clearly faster than r~ d . In terms of the WH model, this 4 He-Vycor 
system falls into the SRC regime (jl.lj) . so theoretical predictions agree with experiments. 
Alternatively, Zassenhaus and Reppy [ [2U] , reported calorimetric studies for this 4 He- 
Vycor system, showing that the singular part of specific heat curves fit well a A-like curve 
with the same exponent a as bulk SF 4 He. Further experiments, using porous gold to 
confine SF 4 He, instead of Vycor, confirm these results [I21|. Microchannels within these 
structures are similar to those of Vycor, but porosities are greater. 

The challenge to theory comes in the case of 4 He confined in silica aerogels (AE) and 
xerogels (XE), where authors observed Cae ~ 0.81 and £ XE ~ 0.89, respectively, larger than 
the bulk exponent. After the WH model, these results suggest that LRC distributions of 
defects must be present within AE and XE. In practice, correlation functions are measured 
by means of (neutron/X-ray) scattered intensities [Ej, thus algebraic (power-law) regimes 
on these plots correspond to algebraic regimes on g(r). At these regimes, objects present 
self-similar (fractal) structure. Thus, provided that the condition (|1.2|) is fulfilled by the 
exponent (a), critical exponents of the system will be modified by disorder. 

Nevertheless, authors argue that silica AE certainly exhibit fractal regimes for several 
length scales, but they point out that beyond a finite cutoff A AE curves enter homogeneous 
SRC regimes. Thus, close enough to the critical point, where the correlation length of SF 
4 He diverges (£ ~ t~ u ), disorder may appear homogeneous at the typical length scale of 
the system, and exponents must cross over to bulk values. For instance, in the case of 95% 
porous aerogels A AE ~ 150nm, and authors estimate £(t) ~ 480nm already at t ~ 10 -4 , 
but although they run experiments up to t — 10~ 5 , the crossover never appearsf^J. They 
actually rule out an apparent violation to Harris criterion, but leave the open question 
about the LRC character of disorder, and propose the existence of LRC within aerogels, 
that cannot be observed through conventional techniques [Cj. 

An attempt to explain these changes on the critical behaviour of SF 4 He will be dis- 
cussed below, in Section |3J showing that well defined LRC structures actually exist within 
aerogels. In the XE case, instead, a clear explanation to changes in critical exponents does 
not yet exist to date. As aerogels, these structures are created through silica sol-gel aggre- 
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gation, but microstructures (strands) are broken and reorganized at the drying stage of 
the process, resulting in more compact (less porous) materials. It has been suggested [122] 
that, close to the flocculation-percolation transition in the aggregation process, strands 
present the same fractal dimension as the LRC structure, but comparative studies do not 
yet reveal self-similarity in xerogels[!23j. 

On the other hand, the liquid-vapor (LV) critical point, known to belong to the 3DIS 
universality class [EI]' gives another scope into the influence of aerogel-like disorder on 
phase transitions. In this case, any type of disorder, correlated or not, is relevant for 
criticality of the pure system. After Harris criterion [[I], a positive specific heat exponent 
for the pure system (a 3DIS ~ 0.109) makes random SRC disorder already relevant. In 
addition, theory predicts [|H[S]> that the O(l) universality class is also subject to changes 
in criticality under a LRC distribution of defects. 

As we quote [ El Ed and show below, both SRC and LRC structures are present 
within aerogels. Thus, the 3DIS confined in aerogel-like distributions of defects (AEIS) 
may be subject to these two competing effects. Results presented by Wong et al. [125]. 
show that the LV critical point of fluid 4 He in aerogels exhibits specific heat curves with 
a finite peak at T c , which suggests that a AE is < 0, clearly different from the pure system 
exponent. In these and further experiments, using N2 instead [I26j. authors also report on 
the order parameter exponent concluding that, within error bars, it is indistinguishable 
from the pure system exponent. Theoretical predictions on the 3DIS with defects, give 
second order corrections to magnetic exponents respect to those of the pure system. In 
addition, the coexistence curve is about 12 times narrower than that of the pure system [ 
I25| . so it is not surprising that changes were so difficult to detect. 

The only trace then left of an effect of aerogel-like disorder on the LV critical point 
is that concerning the finite peak of the specific heat. In a recent paper [EZj> Paredes 
and Vasquez show that this finite peak is consitent with a shift in the correlation length 
exponent to a value close to 2/3, which is closer to the corresponding exponent for the 
randomly diluted Ising system (RDIS) [ [2H1 00] j which is SRC at concentrations of 
disorder c below the percolation critical value p c . However, it is hard to reach definitive 
asymptotic values for the exponents, due to an apparent oscillating approach to the fixed 
point. The same question emerges thus from both, experiments and simulations: why 
criticality would shift to the SRC and not to the LRC fixed point? 

3. Correlations within aerogel-like distributions of disorder: DLCA 

In the sol-gel process of aerogel construction, silica dust is suspended in a solvent 
(sol phase), which allows diffusion and cluster-cluster aggregation to take place. Once the 
process has finished and the gel has been dryed, detailed analyses of structure reveal that 
these objects are neither fractal nor LRC in the WH sense (|1.2|) : only finite regions of 
power-law scaling are observed in scattering intensities from these objects [E]- However, 
at a given time t g of the process, one of the clusters spans all the space of the flask, the 
gelling cluster (GC). Although it has not been tested experimentally, evidence exists that 
these GC are fractal at all scales [E] and LRC in the WH sense [E]- A numerical tool, 
the diffusion limited cluster- cluster aggregation algorithm (DLCA), developed in early 80's 
separately by Meakin[|^| and Kolb, Botet, and Jullien [El], has proved to reproduce well 
the geometrical features of real aerogels [113], 

As silica dust in suspension, in DLCA, particles that initially occupy random posi- 
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tions in a 3D box, are allowed to perform independent Brownian motions, and then kept 
together whenever two of them enter in contact. Undergoing this process, monomers and 
aggregates continue to move randomly, and attach themselves to each other until a unique 
cluster has finally been formed. In the on-lattice version of the algorithm, under periodic 
boundary conditions, cl? particles occupy random sites in a simple cubic L 3 box, being 
c the concentration and <p = 1 — c the porosity. A mass dependent diffusive constant is 
taken into account through a probability p ~ m _<r for an aggregate of mass m to move, 
where ? is a positive mobility exponent. 

As in experiments, the GC is defined as the first cluster to reach opposites sides of 
the box in any direction. Right after the GC is built, many other smaller clusters (is- 
lands) continue to diffuse, and finally attach themselves at random sites to the GC. In the 
fabrication of an aerogel, the resulting structure is preserved by hypercritical evaporation 
of the solvent [ 12], so the GC matrix is kept almost untouched, with randomly dis- 
tributed strands attached to it. Islands represent the SRC subset within the whole cluster. 

Fig. [21 shows the geometrical features of 
different subsets of DLCA clusters, built in 
L = 128 boxes at concentration c = 0.05. 
Following a box-counting procedure [ 13*3*]. 
the scaling of corresponding average densi- 
ties n(r) = m(r)/r d have been depicted as 
function of the box size r. Average over 32 
realizations is taken. Plots marked as BB- 
DLCA and BBGC correspond to backbones 
for the complete DLCA structure and the 
GC, respectively. Plots corresponding to 
complete DLCA clusters (circles) fast enter 
to the homogeneous regime, due to islands 
attached to random sites at gelling clusters 
(squares): the slow power-law decay of GC 
alone suggest that these objects are LRC, 
but when islands have already reached the 
GC to finally form the whole DLCA struc- 
ture, the LRC regime disappears screened (hidden) by the homogeneous distribution of 
islands. Backbones BBDLCA and BBGC also exhibit power-law long regimes, but as we 
argue just below, exponents larger than 3/2 give a > d = 3, which make these subsets 
enter the SRC class, after WH criterion 

The average mass m(r) scales as r f , being dj the fractal dimension of the object. 
Recall that in the box-counting procedure one only takes non-empty boxes to estimate the 
average. This is equivalent to consider n(r) as a one-point correlation function: considering 
only occupied sites as centers for ci-dimensional cubes, the algorithm effectively measures 
the probability for surrounding sites to be occupied, normalized to r d . This method rends 
a decaying exponent (d — df) for n(r). We have argued[EI that the decaying exponent for 
this one-point correlation function is related to that of the impurity-impurity correlation 
function g(r) ~ r~ a as a = 2(d — df). 

An example for this assertion comes from the percolation model, whose density-density 
correlation function scales as g(r) ~ r -( d ~ 2 +v) = r ~ a a t ^he critical occupation fraction 
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Figure 2. Average density n(r) versus box size 
r < L/2 for different subsets of c = 0.05 
DLCA clusters (L = 128). The dashed line 
represents the slope a/2 = 1.5. 
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p = p c . Using Rushbrooke hyperscaling relation one obtains d — 2 + r] = 2(3 /v, being (3 
the critical exponent for the order parameter Poo- Remind that the order parameter in 
the percolation problem is Poo = ra^/N , where moo is the mass (size) of the percolating 
cluster at p c . This quantity scales as P^o ~ L~( d ~ d f\ and being this the order parameter 
for the percolation transition [ 33, its exponent is defined as d — df = (3 /v. In this case, it 
follows that the relation stated above, a = 2(d — df) seems to hold [E2]- 

For the complete DLCA cluster, the crossover to a homogeneous regime occurs at 
r <C L/2, while the fractal regime for backbones seems to span all scales. In the WH 
scheme, however, backbones present the fast decaying SRC regime comparing the 

corresponding n(r) with the slope 1.5 (dashed line). The long-range decay for the GC 
seems slower than r~ d , which enters in the LRC regime (|1,2|) . This case is more complex 
as it becomes homogeneous at r ~ L/2, probably because some strands are already ran- 
domly attached to it at the gelation time t g . To show that GC are in effect fractal at the 
thermodynamical limit L — > oo, and complete DLCA are not, plots of n(r) for different 
lattice sizes L = 32 — 128 are shown in Fig. |3 Arrows indicate crossover lengths which, 
for DLCA clusters (a), reach no more than a few lattice constants, independent of L. On 
the other hand, arrows in Fig. |3] (b) indicate that crossover lengths of GCs scale with the 
lattice size, a trace of fractality of these objects at L — » oo. The complete DLCA cluster 
contains this LRC structure within, with numerous islands randomly attached to it, which 
hide correlations in the box-counting numerical method and make the whole structure ap- 
pear homogeneous. This screening effect is unavoidable in experimental measurement of 
correlations (small-angle scattered intensities), resulting in finite fractal regimes for silica 
aerogels, which tend to disappear as density increases [ II 2j. For simulated aerogels, finite 
size modifies the long-range behaviour of the correlation function. This effect is taken into 
account through the stretched exponential expression proposed by authors in Refs. [l22"l 135] 
n(r) = c g + A r~( d ~ d f^ exp[— {r/£,) s ], where c g is the actual density of the GC (which tends 
to zero as L — > oo) and £ is the characteristic radius of gyration (cuttoff). The parameter 
5 > 1 describes the faster crossover to the homogeneous regime characteristic of finite 
lattices. The dashed curve in Fig. |3] (b) represents the fit to this expression for L = 128, 
giving the estimate df ~ 1.7 (thus a ~ 2.6) for the fractal dimension of the GC, with a 
cutoff at £ 35. 

4. 0(m)-vector models with aerogel-like distribution of defects 

Monte Carlo simulations of 0(m)-vector models consider nearest neighbor interactions 
between spins, placed on a 3D simple cubic lattice of size L (1 < i < L 3 ), with periodical 
boundary conditions. When defects are present, these interactions are described by the 
Hamiltonian: 

^ = -^J>^A'4 (4-1) 

where K = J/kT is the coupling (inverse temperature, with k = J = 1), and = 1 if a 
spin occupies the site, or = otherwise. Depending on the internal dimension (m) of the 
order parameter, spins are taken as (pi = ±1 for the 3DIS (m = 1), and 4>i = (cos 9i, sin 9i) 
for the 3DXY (m = 2) model. For each independent disorder realization, an aerogel is first 
built by the DLCA on-lattice algorithm described above [E2]> taking each site from the 
complete final cluster as a defects. Aerogel pores, i. e., N = (1 — c)L 3 empty sites left in 
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the lattice, are then filled with spins, being c is the concentration of defects. Monte Carlo 
sweeps (MCS) are then taken as follows, using Hamiltonian Q4.1JI : Wolff cluster update 
algorithm is implemented, taking one sweep as eight consecutive cluster flips, to reduce 
autocorrelations caused by critical slowing down. We quote here the procedure followed 
in simulations of the 3DXY model [ I17j. For each disorder realization, thermalization is 
reached after 2 x 10 4 sweeps, while 2 x 10 6 production sweeps are taken at equilibrium for 
further statistical analyses. A suitable number of disorder realizations have been taken, 
depending on lattice sizes, which run from L = 16 (256) to L = 64 (16). The method 
is rather different in the disordered 3DIS[|27], in that shorter time series (~ 10 3 MCS) 
are taken in statistics, increasing the number of realizations at each size (~ 2 x 10 3 ). A 
detailed analysis on self-averaging of thermodynamical quantities, performed for the 3DIS 
case, shows a lack of self-averaging in the magnetization, suceptibility and specific heat, 
and on the other hand, the energy is weakly self-averaged, with an exponent x ~ 2.58. 
In that work [ |2jj , disorder distributions of thermodynamic observables for simulations of 
the 3DIS in aerogels are shown to be more symetric and narrower than for RDIS. Thus, 
good rough estimates could be found with not so much disorder realizations for the aerogel 
case. For this reason, although such analysis has not been performed for results presented 
here for the 3DXY case, an acceptable agreement with experimental results on the SF 
transition of 4 He is obtained [117]. 

4.1. Observables 

The magnetization squared (order parameter) is calculated as, 

M 2 = [|(0)| 2 ] = [|^^| 2 ], (4-2) 

i 

where brackets represent canonical ensemble averages, while square brackets express av- 
erages over disorder realizations, for fixed L and c. The energy per spin is calculated as 
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the canonical ensemble average: 



£=[{E)] = [(H)/N}. 



(4.3) 



Moments for the energy per spin are estimated at the simulation coupling K Q , as £ n = 
[(E n )], as well as moments M n for the magnetization are estimated for the corresponding 
energy bins. Logarithmic derivatives of moments n = 1,2,4 of magnetization are then 
calculated using the corresponding E — M n covariance in the canonical ensemble: 



dln(M r > 



(M n E) - (E){M n ) 



dK (M n ) 
The specific heat is calculated through fluctuations of the energy per spin: 

C = NK- 2 ((E 2 ) - (E) 2 ) . 



(4.4) 



(4.5) 



The helicity modulus [(T^)] has been determined using the Kubo formula [E3E3, 
averaged over the three directions /t = x, y, z, assuming that the model and disorder 
distributions are isotropic. This quantity is proportional to the superfluid density, as was 
derived by Fisher and colleagues [EH], so it is directly related to the order parameter of 
the SF transition. We use this relation to show that the confined system (AEXY) acts in 
d = 3 external dimension, i. e., aerogel pores are effectively three-dimensional [I17j. 

These quantities have been stored from MC simulations at each realization of disorder 
and, after the application of reweighting techniques, extrapolation of thermodynamical 
quantities of interest are made at couplings K « K Q [ [23 OBI ■ 



4.2. Finite size scaling: thermal exponents 

Maxima from reweighted curves for log- 
arithmic derivatives of moments n = 1, 2, 4 
of the magnetization Q4.4J1 . and the specific 
heat (|4.5[) . as functions of the coupling K, 
are averaged over disorder to estimate ther- 
mal critical exponents. 



4.2.1. 0(2)-vector models: 
AEXY 



3DXY and 



According to the finite-size scaling 
(FSS) ansatz [ |IU] maxima of logarithmic 
derivatives scale with lattice sizes L as: 




n \max 



(4.6) 



Respectively, data for the 3DXY and the 
AEXY systems are depicted in Fig. 0J as 
function of lattice sizes L. Both groups of 
data have been vertically displaced (multi- 
plied by a constant) to make corresponding 
plot coincide at L = 16, which allows mak- 
ing a fit on collected data. The FSS power-law fit (|4.6D for bulk 3DXY, using lattice sizes 



Figure 4. FSS of maxima of logarithmic 
derivatives T> n (n = 1,2,4) for bulk 3DXY and 
the AEXY at concentration c = 0.05 of de- 
fects. Corresponding original plots have been 
vertically displaced to coincide at L = 16. 
Dashed lines are fits to the FSS expresion (|4.ti|l 
on collected data. 
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L > 24 gives l/u pure = 1.483(3), which agrees well with previously reported results [HT]. 
and the exponent £ ~ 0.67 reported for bulk A-transition[^j. The same fit to data for the 
AEXY, at concentration c = 0.05, gives 1/V A exy = 1-29(2), clearly different from the bulk 
value. Fits on Fig. 21 (dashed lines) have been vertically displaced to facilitate comparison 
between the slopes. 

The correlation length exponent f AEXY ~ 0.77 obtained for this confined system agrees 
well with the exponent £ ~ 0.76 obtained in experiments [;9j. Both results, numerical and 
experimental, suggest that a LRC distribution of defects must be affecting the critical 
behaviour of the 3DXY universality class, when confined in aerogel-like structures. After 
the extended criterion (|1.2|) , this structure must have a decaying exponent a < 2/ Vpu re ~ 
2.98. Analyses on the DLCA structure resumed above give a ~ 2.6 for gelling clusters 
of DLCA aerogels at c = 0.05, so the WH condition is fulfilled and the distribution GC 
of disorder could be relevant. Nevertheless, there are other correlated structures within 
DLCA aerogels, the backbones, but given that a > 3 for these structures, they do not 
satisfy the WH condition. Had these structures have a decaying exponent smaller than 
d, Fig. [2 shows that this exponent is already greater than the corresponding exponent 
for the GC. The analysis resumed by WH[ |2] when multiple exponents are present on 
LRC strucutures, i. e., when g(r) = Y2i9i r ~ ai i ^ is the term with smallest exponent 
which rules the critical behaviour for the impure system. Indeed, Prudnikov et al. [ [3] 
predict that for a decaying exponent similar to that reported above for the GC (a ~ 
2.6), the correlation length exponent should be z^ LRC ~ 0.73, so both, experimental and 
numerical results roughly confirm this prediction. It is worth recalling that the whole 
DLCA structure itself is not LRC, but it is a mixture of SRC structures with a clearly LRC 
structure, the GC. As follows from WH results, SRC subsets are irrelevant. Thus, evident 
changes observed in critical exponents for this univerality class suggest that only the LRC 
subset rules the critical behaviour of both systems, SF 4 He and 3DXY interacting spins. 

In addition, we have further determined 
critical exponents for the specific heat and 
the helicity modulus as well. Maxima for 
the specific heat, assuming that a < 0, 
scale as: 

[C\ max ~C O0 + AL a /". (4.7) 

These maxima have been depicted in Fig. [5] 
as function of L, normalized to c(L = 
64) = 1, for bulk 3DXY (circles) and 
the confined AEXY system at c = 0.05 
(squares). Non- linear fits of corresponding 
points to the FSS expresion (|4.7|) . using 
data for L > 28 to avoid corrections to 
scaling, are shown in Fig. by correspond- 
ing lines. The method was developed by 
Schultka and Manousakis[ I42j . to determine the specific heat exponent for the pure 
3DXY system. The corresponding best fits give a/u ~ —0.015 for the pure system, 
and ctjv ~ —0.38 for the confined system. These results agree well with our results 
on correlation length exponents, shown above, after the Josephson hyperscaling rela- 
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Figure 5. FSS of maxima of the specific heat for 
bulk 3DXY (circles) and confined AEXY at con- 
centration c = 0.05 of defects (squares). Lines are 
fits to the FSS expresion (|4.7|l . 
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tion [ I18j in the form 2/v — a/u = d* , which gives d* ~ 3.0 for both systems. Us- 
ing corresponding results on the correlation length exponents, we obtain a ~ —0.010 
for bulk 3DXY, which agrees well with calorimetric studies on SF 4 He[ 1£L The result 
a = —0.29 for the AEXY is lower than the reported experimental value a — 0.57 for 
4 He confined in aerogel pores at the same volume fraction c = 0.05 [^J. However, in 
numerical simulations reported here, the question posed by authors about the apparent 
violation of hyperscaling seems not to emerge from our results on exponents v and a. 
The critical coupling K c has been determined using the crossing of reweighted plots of the 
Binder fourth cumulant for the mag netization f7 4 = 1 - (M 4 )/3(M 2 ) 2 at different lattice 
sizes L, which is universal at K c (independent on L)[l43j. This analysis yields K c ss 0.45416 
for the pure system, while the crossing of average reweighted curves in the AEXY case 
gives K c ps 0.46495. 

Using these values, we estimate the he- 



licity modulus at K c , 
der realizations: 



averaged over disor- 
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[<T A >] ~ L~ v l v . 



(4.8) 
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Fig.Elshows corresponding data for the he- 
licity modulus in the pure 3DXY case (cir- 
cles), and the AEXY case (squares). Lines 
represent fits to the FSS expression 14.81 
The exponent v/v = 1 confirms that the 
effective external dimension for the system 
is d = 3, through the hyper universality re- 
lation v = (d — 2)v [ 39]. This result also 
confirms that the correlation length expo- 
nent we have estimated here corresponds 
to the exponent £ for the SF transition of 
4 He in aerogels. 

Results reported here seem to answer the question quoted above from Refs. [EH Ej: 
numerical evidence has been presented here (Section |3J) that a physically well defined 
hidden LRC distribution of defects exists in aerogels, which explain changes on the critical 
behaviour of SF 4 He when confined in those structures. 



Figure 6. FSS plots of the helicity modulus at 
K c . Bulk 3DXY (circles) and confined AEXY 
at concentration c = 0.05 of defects (squares). 



4.2.2. 0(l)-vector models: 3DIS and AEIS 

We turn now to a brief discussion on changes of the correlation length exponent for 
the 3DIS, when confined in DLCA aerogels at concentration c = 0.2 of defects. A detailed 
analysis is reported elsewhere [ I27j. The method, as stated above, is quite different from 
that used in the AEXY case in that shorter MC time-series are used for extrapolation, 
but total numbers of realizations have been incremented in two orders of magnitude. 
Shorter CPU times, in this case, allows also making simulations at greater lattice sizes, 
being L = 8 — 96 in this case. After extrapolation by reweighting, curves have been 
averaged over disorder for each K, so a unique curve is obtained for each set of data. 
Position of maxima {pseudocritical couplings K*{L)) for the magnetic susceptibility (x), 
and logarithmic derivatives T> n of moments of the magnetization (n = 1,2,4) have been 
estimated through the averaged curves. According to the finite size scaling theory [ 143) . 
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maxima for a given observable O scale as 

K*(0,L) = K c + a{0)L- l l y + c. s. (4.9) 

being K c the critical coupling at the thermodynamical limit, L — ►oo. Coefficients a{0) are 
non-universal, i. e. , dependent on the observable O and details of the system. Corrections 
to scaling (c. s. ) are avoided in the analysis reported below, taking data for large enough 
lattice sizes (L > 40). Points depicted on the inset of Fig. have been calculated through 
local differences K*(Oi,L) — K*(Oj,L) between data from two different observables Oi 
and Oj. This transformation of data allows eliminating the non-singular term (K c ) from 
Eq. Q4.9j) . common by hypothesis to all observables, allowing an easier estimation of the 

Power-law fits have been made to these 
transformed data, giving a rough estimate 
of the (common) FSS exponent 1/V ~ 1.41. 
This approximate has been used to rescale 
the horizontal axis to L -1 ^, as indicated 
in Fig. As expected, a linear behaviour 
is observed for all data, and linear fits to 
the expression (|4.9|) gives just one value 
K c = 0.25857(1) for all plots. A more de- 
tailed analysis on the scaling of logarith- 
mic derivatives at the critical point reveals, 
however, that the correlation length expo- 
nent tends to a value closer to 3/2 [ [22]. The 
correlation length exponent vaeis ^ 0.67 
obtained is clearly greater than that for 
the pure 3DIS and the LV critical point 
{y w 0.63). These results point to a neg- 
ative or zero value for the specific heat ex- 
ponent a. Experiments on the LV critical 
point of 4 He in aerogels [ 125} report specific heat curves that seem to have a finite peak 
at T c , but authors do not estimate the critical exponent a which may either be negative 
or zero as well. After Josephson hyperscaling, this is consistent with our rough result 
1/v < 1.5, tending to values closer to that predicted for the SRC Ising universality class 
than that predicted for the LRC fixed point [IH]. 

We must then conclude, in this case, that the FSS method applied here may not 
be precise enough to determine really asymptotic values for critical exponents, but a 
remarkable difference exist between results approached here and exponents for the pure 
system, as expected after Harris criterion. 

5. Concluding remarks 

We have reported here on the effects of a specific type of disordered distributions of de- 
fects on the critical behaviour of 0(m)-vector models, using the 0(1) and the 0(2) models, 
represented by the 3DIS and the 3DXY models, respectively. A detailed scaling analysis 
has been made to structures created through the DLCA on-lattice algorithm, which repro- 
duces well the geometrical tasks of real silica aerogels. This analysis shows that the whole 



exponent on the leading term aL l l v . 



0.2588 



0.2586 

K 

0.2584- 



0.2582 



0.2580 



0.2578 



Y — 
-O 



1/V = 1.41 



0.002 



-Q— 



— £3- 









— EL- _ 



































X 






□ 


D i 






o 


D 2 






A 


D 4 



£~ 1/v 0.004 



-H— 



0.006 
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positions of maxima for the susceptibility (x), 
and logarithmic derivatives (£>„ ) . The horizon- 
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structure is not, strictly speaking, a fractal: the power-law decaying regime is preserved 
only through a few lattice constants, and even more, it does not scale with lattice sizes. 
A more detailed dissection, shows the existence of a physically well defined subsets within 
the whole strucuture, gelling clusters, that keep the algebraic decaying regime through 
lengths comparable to the lattice size, with a cutoff that scales with the system size. Un- 
der the light of WH theory, the existence of this fractal (LRC) subset already explains 
changes on critical exponents of the 3DXY model, when confined in these aerogel-like sim- 
ulated structures, respect to those of the pure system. Results for the confined (AEXY) 
system roughly confirm those observed for the SF transition of 4 He inside aerogels, and 
explain why a crossover to bulk values is not observed in critical exponents reported from 
experiments, as T — > T c , when the correlation length divergesf^J. Even more intriguing 
are results from the 3DIS confined in the same type of disorder distribution, as exponents 
seem to be affected by the LRC fixed point to finally tend values far from those predicted 
by theory[[Sj, and closer to those of the SRC fixed point. The question about the critical 
behaviour of the AEIS rests open as conclusive results on asymptotic values for critical ex- 
ponents have not yet been obtained. Previous renormalization group flows to fixed points, 
sketched by Weinrib and Halperin[|2], and analyses made by Prudnikov et al. [Ej, point 
out that the approach to the LRC fixed point may be oscillating, which can make very 
dificult to enter a really asymptotic regime applying conventional finite size scaling. 
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Questions and answers 

Q (Rafael Rangel): What are the theoretical predictions for the (possible) 
change of exponents by disorder? 

A : Harris criterion predicts that a weak uncorrelated disorder is not relevant if the 
correlation length exponent v is larger than 2/d for the pure system, or after Josephson 
hyperscaling, a = 2 — dv < 0. For 3DIS case, a > 0, and then short-range-correlated 
disorder could be relevant. Weinrib and Halperin extend this criterion to cases where 
there exist long-range correlations (LRC) in the distribution of defects, as expressed by 
Eqs. 1)1.1(1 and (|1.2j) . It is this extension that allows explaining changes on the critical 
behaviour of the 3DXY model when confined in aerogel-like distributions of defects, and 
changes on critical exponents of the SF transition of 4 He in aerogels as well. 

Q (Dragi Karevski): How does your FSS technique work if there are different 
correlations? 

A : Weinrib and Halperin [|2j already considered the case where there exist terms with 
different decaying exponents in the correlation function g(r), as was quoted in Subsection 
14.2.11 the critical behaviour couples to the term with the lowest exponent. Indeed, in 
the model disorder we are considering in simulations presented here, DLCA, there exist 
several subsets of the whole cluster which present different decaying regimes. Results point 
to state that, in the 3DXY confined case (AEXY), it is the gelling cluster (GC) that rules 
the critical behaviour, as it is the LRC subset with the lowest exponent (see Fig. |2J). 
The AEIS case is by far more interesting in that, competing effects are present, as SRC 
disorder is also relevant. Being the LRC fixed point marginal, it seems that it is the SRC 
fixed point that finally governs the critical behaviour, i. e., the effective exponent 1/v turns 
from being attracted to 1/v ~ 1.4 (predicted by Prudnikov et al. for the LRC fixed 
point) to 1/v < 3/2, closer to the SRC fixed point value. Nevertheless, even if only one 
model LRC distribution of defects were present, the application of FSS techniques must 
consider that the convergence of the correlation function to the final decaying regime 
j^-(d-d f ) w jvj depend on the lattice size L, so the convergence of effective exponents to 
critical values will be ruled by an exponent that asymptotically tends to a stable value. 
This task of the correlation function may cause succesive crossover regimes which appear 
in the convergence of effective exponents to asymptotic critical exponents. 

Q (Dragi Karevski): How does the FSS with V = L d (and not with L) hold? 

A : The FSS theory is based on the homogeneity of the free energy density, which 
allows writing it in terms of the scaling variable L/£ = Lt u , which at its time may be 
redimensioned to tL~^l v . The external dimension d of the system appears only in corre- 
sponding hyperscaling relations, such as those of Josephson 2/v — a/v = d and Rushbrooke 
7/1/ + 2/9/i/ = d. 

Q (Wolfhard Janke): What is the dimension dt of your fractal? 

A : The whole DLCA clusters in our model disorder are not fractal, as they enter ho- 
mogeneous regime at very few lattice constants, but the well defined GC have correlations 
that seem to span all scales at the thermodynamical limit. The best estimate we have of 
the fractal dimension for these objects is df ~ 1.7. 

Q (Yurii Holovatch): How is c related to a? 

A : As the concentration of defects increases, the fractal dimension of the GC increases 
as well, as stated in Ref. [H3j. so the exponent a decreases. 
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